Granular discharge and clogging for tilted hoppers 
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We measure the flux of spherical glass beads through a hole as a systematic function of both tilt 
angle and hole diameter, for two different size beads. The discharge increases with hole diameter 
in accord with the Beverloo relation for both horizontal and vertical holes, but in the latter case 
with a larger small-hole cutoff. For large holes the flux decreases linearly in cosine of the tilt angle, 
vanishing smoothly somewhat below the angle of repose. For small holes it vanishes abruptly at a 
smaller angle. The conditions for zero flux are discussed in the context of a clogging phase diagram 
of flow state vs tilt angle and ratio of hole to grain size. 

PACS numbers: 45.70.-n 



I. INTRODUCTION 

The flow of granular materials is of widespread practi- 
cal and fundamental interest. One challenge 
to understanding and controlling behavior is that the re- 
sponse is nonlinear, with a forcing threshold below which 
the medium is static. Furthermore, just above threshold 
the response may be intermittent even though the forc- 
ing is steady. Familiar examples include avalanches down 
the surface of a heap as well as gravity-driven discharge 
from a horizontal hole at the bottom of a deep container 
or "silo" . For the latter, the mass discharged per unit 
time is given by the "Beverloo" relation: 

W ^Cpb9^'\D-kdf/\ (1) 

where pb is the density of the bulk granular medium, 
g — 980 cm/s^, D is the hole diameter, d is the grain 
diameter, and C and k are dimensionless fitting parame- 
ters 0. By contrast with viscous fluids, the discharge of 
grains is independent of filling depth. This may be un- 
derstood in terms of the classic Janssen argument that 
pressure vs depth approaches a constant for a deep con- 
tainer due to support of the weight of the medium by 
frictional contacts with the sidewalls. While the Beverloo 
relation is supported by a large body of work, as reviewed 
by Nedderman and Savage et al. discrepancies of up 
to forty percent have been recently reported when the 
hole size is increased more widely than usual Q . Typical 
ranges for the numerical constants are 0.5 < C < 0.7 and 
1.2 < fc < 3, depending on grain shape and friction. 

The Beverloo equation implies the existence of a 
threshold hole diameter, kd, of a few grains across, be- 
low which the flux vanishes. Just above this threshold, 
the flow is subject to intermittent clogging [ol. fiol. [Til. [T3| . 
Even far above threshold, the response may not be steady 
in that the Beverloo form is often interpreted in terms 
of intermittent formation and breakup of arches across 
the hole. In particular, grains in a freshly-broken arch 
free-fall through a distance set by hole size and hence 
emerge with speed v {gDy^^ and mass per unit time 
W -iPh^y^-D^/^^ Density waves [H, and ticking 
[IE Ei Tl7| are also examples of unsteady response in 
gravity-driven discharge, but where air plays a role. 




FIG. 1: Schematic cross-section of the discharge apparatus 
tilted at various angles, as labeled. The container wall is in- 
dicated by the heavy black lines, with a circular hole through 
which the grains emerge. The granular material is indicated 
by the gray shaded region. Below some boundary indicated 
by the thin dashed line the grains may be stagnant. Transient 
arches may extend across the hole, as indicated by the thick 
dark gray curves. 



To develop a deeper microscopic understanding of 
granular discharge, it seems important to grapple with 
the unsteadiness represented by intermittent clogging 
and arch formation / breakup. One experimental ap- 
proach is to vibrate the system, in order to fluidize and 
break arches as well as to introduce a competing time 
scale [H, llE HO, HH In this paper, our approach is 
to interfere with the usual arch formation and breakup 
by tilting the container as depicted in Fig. [l] so that the 
plane of the hole is inclined by an angle 6 away from 
horizontal. According to the "free-fall arch" picture, one 
might expect the discharge rate to decrease according to 
the reduced horizontal projection of the hole, cos6', and 
to vanish at a tilt angle less than ninety degrees where 
the projection falls below a nonzero threshold. 

There are few prior experiments on granular flux from 
holes that are not horizontal, as noted in an article con- 
cernin g th e use of vertical slots for constructing a flow 
meter [23|. Ref. [l^l predates Beverloo and reports flow 
rates for two media and three hole diameters at inclina- 
tion angles of {0, 30, 60, 90} degrees. The results are 
claimed to be linear in the cosine of the tilt angle; how- 
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ever, we concur with statements in Ref. 23| that the data 
are too sparse and uncertain to demonstrate the proposed 
form. The definitive review by Nedderman et al. Q cites 
Ref. [25II as a "prehminary investigation which comes to 
no clear conclusion" regarding discharge through a ver- 
tical hole; it also cites Ref. [241 but only regarding hori- 
zontal holes. Ref. [2^ reports that the discharge rate for 
a vertical hole at the very bottom of a sidewall scales as 
diameter to a power between 2.5 and 2.8, and that the 
ratio of vertical to horizontal discharge rates is between 
0.37 and 0.50. Ref. [13] reports on velocity fields, but 
not discharge rates, for conical hoppers tilted up to 23°. 
Thus, the observations reported here concern a relatively 
unexplored effect. 



II. MATERIALS AND METHODS 

The granular medium consists of spherical glass beads, 
with two different diameters: d — 0.30 ± 0.05 mm and 
d — 0.9 ± 0.1 mm. Both have bulk density pi, = 
1.53 ± 0.01 g/cui^ and draining angle of repose 9r = 24°. 
Two different types of container are used. The first type 
is a steel can, with 10 cm diameter, 12 cm height, and 
0.25 mm wall thickness. Holes are drilled in three differ- 
ent locations: in the bottom at center, in the bottom at 
2 cm from the side, and in the side at 2 cm above the 
bottom. The second container type is square Aluminum 
tubing, with 7x7 cm^ inner cross section, 31 cm height, 
and 3 mm wall thickness. One to four holes are drilled in 
each side in a staggered arrangement near the bottom, no 
closer than 2 cm from any edge or from each other, and 
are countersunk with a 120° chamfer. The containers are 
open at the top, so there is no back-flow of air into the 
hole to replace the loss of granular material. The contain- 
ers are grounded to prevent electrostatic charging, and 
are mounted by a chain clamp with changeable orienta- 
tion on a ring stand. The tilt angle 9 of the plane of the 
hole away from horizontal is measured with a plumb bob 
and protractor; 6* = 0° corresponds to a horizontal hole, 
9 — 90° corresponds to a vertical hole, and 9 > 90° corre- 
sponds to discharge with an upward angle of 9 — 90°; see 
Fig. [T] for a sketch of the system at four different tilt an- 
gles. Flow is impossible for 9 > 180° — 9r = 156°, where 
the granular medium loses contact with the boundary 
around the hole. 

Discharge rates are measured by weighing the material 
collected during a timed interval ranging from several sec- 
onds for fast flows, to several minutes for slow flows, as 
follows. First the hole is covered with a piece of paper 
and the granular medium is poured from a beaker into 
the container, usually to its full height. Then the cov- 
ering is removed and flow allowed to proceed for a sev- 
eral seconds or more. Next a beaker is inserted into the 
discharge stream while simultaneously starting a timer. 
Finally, the beaker is removed from the stream while si- 
multaneously stopping the timer. Statistical uncertainty 
in discharge rate is typically one to ten percent, as re- 
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FIG. 2: (Color online) Discharge rate vs hole diameter for (a) 
horizontal and (b) vertical holes; the ratio is shown in (c). 
Red circles (blue diamonds) are for d — 0.3 mm (0.9 mm) 
diameter glass beads. Open symbols are for containers made 
from AI tubes of square cross section; solid symbols are for 
steel cans. Solid curves are two-parameter fits to the Beverloo 
form, as specified; the ratio of these fits is shown in (c). 



fleeted by the size of the scatter in the data for runs taken 
under identical or similar conditions; error bars are not 
displayed since they are smaller or comparable to symbol 
size. 



III. DISCHARGE RATE 

We begin by considering the masses per unit time, Wq 
and 14^90 , respectively discharged through horizontal and 
vertical holes. Initial attention is restricted to hole sizes 
that are large enough not to exhibit clogging or to require 
tapping to initiate flow. The first observation is that the 
discharge rates do not depend on the filling depth, which 
was varied from the full container height to about one- 
third of the container height for both the steel cans and 
the Aluminum tubes. Care is thus taken to ensure that 
measurements are not corrupted by any change in dis- 
charge rate that may arise as the medium fully empties 
from the container and the upper free surface comes near 
the orifice. This well-know behavior contrasts with liquid 
discharge, and is explained by the classic Jansscn argu- 
ment. Furthermore, the discharge rates do not depend 
on details of the container geometry, which was varied in 
terms of container shape and hole location as described 
above. This is also evident in Figs. [Ha,b), where dis- 
charge rates for the two container types are plotted ver- 
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sus hole diameter, for both horizontal and vertical holes 
and for both grain sizes. 

The observed discharge rates all increase with hole 
size, more rapidly for smaller holes, and appear to ap- 
proach a power-law for larger holes. Fits of discharge 
data to the Beverloo relation are included in Figs. [21Ja,b) 
as solid curves. The agreement is good and furthermore 
the numerical constants are independent of bead size: 
{Co = 0.573 ± 0.003, ko = 1.6 ± 0.1} for horizontal holes 
and {Cgo = 0.269 ± 0.005, fcgo = 2.3 ± 0.2} for verti- 
cal holes. Note that, therefore, the two curves in each 
plot represent one simultaneous fit. The fitting parame- 
ters for horizontal holes in Fig. [Ifa) are in accord with 
prior work 7] . The fitting parameters for vertical holes in 
Fig. ^h) indicate both that the zero-flux threshold hole 
diameter is larger than for horizontal holes, kg^d > kgd, 
and that for very large holes the ratio of vertical to hor- 
izontal discharge rates becomes constant. This can be 
seen directly in the discharge ratio, Wgo/VFo, plotted in 
Fig. IHc). The observed ratios increase with hole size 
and fall in the range 0.3 — 0.5 reported in Ref. [26] • Note, 
however, that the asymptotic value Cqq/Cq — 0.47 from 
the Beverloo fits is not attained even for the largest holes 
studied. 

We have no theoretical explanation for the observed 
diameter-dependence of the discharge rate for vertical 
holes, other than that the basic scale must be set dimen- 
sionally as pbg^^^D^^^ with a numerical prefactor that 
depends on inclination. The fact that the Beverloo form 
successfully describes both horizontal and vertical dis- 
charge rates raises questions about the physical inter- 
pretation in terms of free-fall through a distance set by 
hole size, since the component of grain velocity that con- 
tributes to discharge is perpendicular to gravity for verti- 
cal holes. Transient arches could extend from the top of 
the hole to the solid packing of stagnant grains below a 
boundary, perhaps set by the angle of repose, as sketched 
in Fig. [TJ however, upon breaking, the grains in such an 
arch cannot escape at free-fall speed. Any free-fall will be 
terminated by collision with either the stagnant pile or 
with other grains in a flow field that curves toward the 
exit. It would be interesting to map out the flow field 
and the time-averaged density for tilted hoppers, even in 
a quasi- two dimensional system as in Refs. [l9l. [20l. I2R . 
to directly investigate such behavior. In addition, a fur- 
ther question is raised for how to interpret the different 
zero-flux threshold hole diameters, kod 7^ fcgod, in terms 
of an "empty annulus" where grain centers may not pass 
Q because that concept is independent of hole orienta- 
tion. 

Since there is no theory to test, we now examine trends 
empirically two different ways. First, flux data are ob- 
tained for a wide range of tilt angles, not just zero and 
ninety degrees, for both size grains and for three to six 
different hole diameters. We find that the Beverloo equa- 
tion gives satisfactory fits to flux vs diameter in all cases. 
Results for the fitting parameters, C and k, are plotted vs 
cosine of the tilt angle in Fig. [31 The values displayed are 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 




- 




(a) _ 


- 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 

1 
1 


1 1 1 1 
(b) _ 








1 1 \ ' 

, . . , 1 . , , . 1 , . , , 




- 


1 1 


1 1 1^ 






(c) " 









O^t^fi — , — , — I — , — , — , — , — I — , — , — , — , — \ — , — , — , — , — I 
-1 -0.5 0.5 1 



cose 

FIG. 3; (Color online) Dimensionless parameters, C and fe, 
obtained by fitting flux vs diameter data to the Beverloo re- 
lation, W = C pbg^^^{D — kd)^''^, for holes oriented at angle 6 
away from horizontal. Solid symbols are for the fits shown in 
Fig. [2] and open symbols are from fits to data shown in Figs.|4]- 
[S] Error bars are based on both uncertainty in fits as well as 
on the difference of results for the two grain sizes. The ver- 
tical dashed line and shaded region indicate cos 6 < ~ cos Or 
where discharge is impossible. The solid lines in (a) and (c) 
are the best fits to a line; the solid curve in (b) is the best fit 
to a parabola with minimum at cos 6 = 1. 



an average for both d = 300 /im and d = 900 /im diam- 
eter grains. Note in the top plot that C decreases with 
tilt angle and appears to be a linear function of cos6'. 
Furthermore, it extrapolates to zero at — cos^^, below 
which no flow is possible. Note in the middle plot that 
k increases with tilt angle, though the functional form 
is not as clear. It appears to depart quadratically from 
the zero-angle value. And since C apparently vanishes at 
Or, it is natural to speculate that k diverges at the same 
angle. This possibility is reinforced in the bottom plot, 
in which 1/fc^ appears to be a linear function of cos 9 and 
to vanish at — cos Or ■ 

Alternatively, the tilt-dependence of the flux can be 
considered empirically without using the Beverloo rela- 
tion. Speciflcally, for a given hole size, we examine how 
tilting causes the flux to decrease from a maximum at 
zero-angle, Wq. The simplest hypothesis would be re- 
duction according to the projected horizontal area of the 
hole, i.e. W = Wq cos 9. So we plot discharge rates for 
various hole diameters vs cosO in Fig. |3]for d = 0.3 mm 
diameter grains. Raw data are shown on a logarithmic 
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FIG. 4: (Color online) (a) Discharge rate vs cosine of tilt 
angle for d = 0.3 mm diameter beads through various hole 
diameters D, as labeled; values of D/d are given in parenthe- 
sis, (b) Discharge rate normalized by value at zero tilt angle, 
with same symbol codes. The vertical dashed line and shaded 
region indicate cos 61 < — cos^r where discharge is impossible. 
The data points with stem lines indicate angles at which the 
flow is subject to clogging. In (a) multiple points at a given 
angle represent data taken with different container types and 
hole locations; in (b) each point represents the average over 
different container/hole geometries, and the black line repre- 
sents the best linear fit to the large-hole data. 



scale in part (a), while the ratio W/Wq of flux to the value 
at = is shown on a linear scale in part (b) . There are 
several interesting features in these plots. First, as shown 
by multiple data points at a given angle, the results are 
independent of container geometry: square tubes and 
cylindrical cans with different hole placements all have 
the same discharge rates. Second, the data sets appear 
nearly parallel on the logarithmic scale of part (a), and 
hence are nearly proportional. However, the data sets are 
not truly proportional to one another since good collapse 
is not found in part (b). For large holes, though, the data 
appear to approach a common linear dependence on co- 
sine of tilt angle: W/Wq = (cos 6* -I- a)/(l-|- a) where the 
fitting parameter is a = 0.78 ± 0.01. We stress that this 
form is linear in, but not proportional to, cos 9; therefore, 
the flux is not proportional to the horizontal projected 
area. Moreover, there is a nonzero flux even for 9 > 90°, 
when the hole is tilted past vertical and the unit normal 
vector to the hole has an upwards component. Naturally 
no flow is possible for 9 > 180° — 9r = 156°, when the 
medium loses contact with the boundary surrounding the 
hole. This argument underlies the speculation in Ref. [l^l 
that the value of a should be cos^r- In fact, the flux van- 
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FIG. 5: (Color online) (a) Discharge rate vs cosine of tilt 
angle for d — 0.9 mm diameter beads through various hole 
diameters D, as labeled, in the square Al tube containers; 
values of D/d are given in parenthesis, (b) Discharge rate nor- 
malized by value at zero tilt angle, with same symbol codes. 
The vertical dashed line and shaded region indicate the re- 
gion cos < — cos 9r where discharge is impossible. The data 
points with stem lines indicate angles at which the flow is 
subject to clogging. The black line in (b) represents the best 
linear fit to the large-hole data. 



ishes below this bound at an angle arccos(— 0.78) = 140°. 
All these features, including the value of a, can also be 
seen in Fig. [5] for d = 0.9 mm grains. Theoretical guid- 
ance would be helpful in determining whether the best 
description of discharge rates should be formulated in 
terms of Bcvcrloo fit parameters, as in Fig. [31 or in terms 
of reduction from zero-tilt discharge, as in Figs.[l][51 

We note that the Beverloo parameters in Fig. [3] were 
obtained by fits to the raw data displayed in the top plots 
of Figs. [Hand [5l reorganized as a function of diameter at 
a given angle. 



IV. CLOGGING 

In this final section we consider whether or not the 
system is clogged or fiowing, and how the state of be- 
havior depends on hole size and tilt angle. One might 
suppose that the zero-fiux threshold hole diameter, kd, 
found from fits of the Beverloo form to flux data at dif- 
ferent tilt angles would demarcate the phase boundary 
between clogged and flowing; however, this is not the 
case. To investigate, we first note that some features 
of clogging are already displayed in the plots of W and 
W/Wq vs tilt angle in Figs. Specifically, stem lines 
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FIG. 6: (Color online) Clogging phase diagram of discharge 
state vs tilt angle and ratio of hole to grain diameter. Circles 
are for d — 0.3 mm grains; triangles are for d = 0.9 mm 
grains; open symbols are for stop angles at which flow ceases 
as the tilt angle is increased; solid symbols are for start angles 
at which flow commences as the tilt angle is decreased. The 
horizontal dashed line and shaded region indicate the region 
9 > 180° — Or = 156° where discharge is impossible. The 
solid black curves are guides to the eye that rise abruptly 
from zero, bracket the transition, and approach 9r for large 
holes. The dotted green curve represents the normalized zero- 
flux hole diameter, k, from fits to the Beverloo relation, shown 
previously in Fig. |3jb) . 



extending from the data points down to zero are used 
to indicate conditions where the system was observed to 
clog. In these cases, smooth steady discharge proceeds 
for some extended time interval long enough to measure 
flux, but then suddenly and unpredictably stops. With 
gentle tapping or poking, the clog can be broken and an 
interval of steady flow can be restarted. Visual inspec- 
tion reveals that clogging is not caused by impurities or 
larger grains in the tail of the size distribution that block 
the hole. As demonstrated in Figs. HJH we find no such 
intermittent clogging for sufficiently large holes; instead, 
with increasing tilt angle, the discharge rate decreases 
continuously toward zero. For smaller holes, intermit- 
tent clogging happens at nonzero flux over a range of tilt 
angles. Thus the transition between zero and nonzero 
steady flux is discontinuous, but not necessarily sharp. 

As a first attempt to locate the transition and ascertain 
its sharpness we perform the following experiments. For 
a given hole size, first we begin at a small tilt angle in the 
fiowing regime and then slowly increase the tilt angle un- 
til a clog first occurs; we repeat 2-5 times and record the 
average "stop" angle. Next we begin at a large tilt angle 
in the clogged regime and then slowly decrease the tilt 
angle until until flow commences; we repeat 2-5 times and 
record the average "start" angle. Tilting is done by hand 
at a constant rate of roughly one degree per second; this 
is slow enough, and the motion is smooth enough, that 
acceleration and jerk do not influence behavior. Given 
these steady tilt rates, and the absence of a dwell time 
at any particular angle, it takes less than a minute to go 
from 0° to stop or from 180° to start. Note that inter- 
mittent flow with a lifetime of several minutes or more 



will appear continuous. While the stop and start angles 
may thus depend on tilt rate, the average ought to be less 
sensitive; this concurs with informal trials at slightly dif- 
ferent speeds and with actual data where the tilt rate is 
only imperfectly controlled by hand. Thus we construct 
a "clogging" phase diagram in Fig. [6] by plotting the stop 
and start angle data as a function of dimensionless hole 
diameter D j d. For very small holes where no flow occurs 
without tapping, we plot the stop and start angles at 
zero. For holes larger than about D/d = A the transition 
angles rise abruptly from zero. They increase monotoni- 
cally and appear to approach 180° — 9r = 156° for very 
large holes. Note that the stop angle data for different 
size grains exhibit good collapse for all hole diameters. 
By contrast the start angle data appear to collapse well 
only for large holes. 

The sharpness of the transition between clogged and 
flowing states may be gauged from Fig. [S] by the differ- 
ence in start and stop angles. The spread in these angles 
is greatest for small holes, just above D/d ^ 4, and de- 
creases steadily for larger holes. We speculate that the 
increasing sharpness of the transition may be connected 
with the simulation result that the distribution of pack- 
ing densities at which a granular system jams becomes 
narrower as the system size increases js^l • Here, for clog- 
ging to occur, the sample need not be jammed every- 
where - only in a volume Vhoio ~ over the outlet. 
Flow proceeds steadily only until a grain conflguration 
arises in T4ioie that is jammed. For smaller holes, more 
packing fractions exist that can jam and the system is 
more susceptible to clogging. For large enough holes, no 
such configurations exists and the system flows freely or 
not at all. It would be interesting to explore this possi- 
ble connection to Ref. [30| with real-time measurements 
of packing density in I4ioic- It would also be interesting 
to map out the clogged/flowing phase boundary by al- 
ternative means, such as contour plots of average flow 
duration and mass discharged (sij . 

Lastly we note that the shape of the clogged/flowing 
phase boundary may be described by empirical curves of 
the form 9 = 156° (1 - Ad/D)" for D > Ad. The up- 
per solid curve in Fig. [5] through the start angles corre- 
sponds to A = 4.0 and a = 1/4; the lower curve through 
the stop angles corresponds to A = 4.8 and a — 1/3. 
These approach 180 — 0^ for very large holes; however, 
we caution that the actual data do not reach this limit 
and could well saturate at a smaller angle. Inversion gives 
the hole diameter for crossover from clogged to flowing as 
D = Ad/[l-(6'/156°)i/"]. This is quite different from the 
zero-flux threshold hole diameter kd deduced from fits to 
the Beverloo equation, shown previously in Fig. [Sjs and 
now plotted on the clogging phase diagram as a dashed 
curve. Comparison shows that that the clogged/flowing 
hole diameter is about three times larger than the Bever- 
loo zero-flux threshold hole diameter. Thus the suscep- 
tibility to clogging does not directly correspond to the 
Beverloo equation and the vanishing of flux. It would be 
interesting to obtain discharge data for holes below clog- 
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ging, as the flow duration and flux both approach zero, to 
see if deviation from the Beverloo form can be detected. 



V. CONCLUSION 

In this paper we reported on discharge rates and clog- 
ging behavior for glass beads and circular apertures as a 
function of both hole size and inclination angle. Exten- 
sive and syHtcrnatic variation of the latter serve to fill a 
particularly unexplored void in the literature. Our dis- 
charge results shed new light on the Beverloo relation, 
particularly the free-fall arch and empty-annulus inter- 
pretations, as well as its validity for small holes. Our 
clogging results emphasize the need for theoretical con- 



sideration of fluctuation and jamming effects, especially 
for slow flows. Altogether, our experiments round out 
the phenomenology of granular discharge, highlight the 
unusual and elusive mechanics of granular materials, and 
suggest specific further lines of research. 

Acknowledgments 

We thank T. Brzinski for experimental assistance, and 
S. R. Nagel for an inspiring series of experiments and 
simulations on granular media and jamming. Our work 
was supported by the NSF through grant DMR-0704147 
and by the University of Pennsylvania through its work- 
study financial aid program for undergraduate students. 



[1] R. M. Noddcrman, "Statics and kinematics of granular 
materials" (Cambridge University, NY, 1992). 

[2] F. J. Muzzio, T. Shinbrot, and B. J. Glasser, Powder 
Teclinology 124, 1 (2002). 

[3] H. M. Jaeger, S. R. Nagel, and R. R Behringer, Rev. 
Mod. Rhys. 68, 1259 (1996). 

[4] J. Duran, "Sands, powders, & grains: An introduction to 
the physics of granular materials" (Springer, NY, 2000). 

[5] "Jamming and Rheology: Constrained Dynamics on Mi- 
croscopic and Macroscopic Scales" edited by A. J. Liu 
and S. R. Nagel, (Taylor & Francis, NY 2001). 

[6] W. A. Beverloo, H. A. Leniger, and J. Van de Velde, 
Chem. Eng. Sci. 15, 260 (1961). 

[7] R. M. Nedderman, U. Tiiziin, S. B. Savage, and G. T. 
Houlsby, Chem. Eng. Sci. 37, 1597 (1982). 

[8] C. Mankoc, A. Janda, R. Arevalo, J.M. Pastor, I. 
Zuriguel, A. Garcimartin, and D. Maza, Gran. Matt. 9, 
407 (2007). 

[9] S.S. Manna and H.J. Herrmann, Eur. Rhys. J. E 1, 341 

(2000) 

[10] I. Zuriguel, L. A. Pugnaloni, A. Garcimartin, and D. 
Maza, Rhys. Rev. E 68, 030301R (2003). 

[11] I. Zuriguel, A. Garcimartin, D. Maza, L. A. Pugnaloni, 
and J. M. Pastor, Rhys. Rev. E 71, 051303 (2005). 

[12] S. Saraf and S. V. Franklin, Bull. Am. Rhys. Soc. 54, 
J14.00005 (2009). 

[13] G.W. Baxter, R.P. Behringer, T. Fagert, and G.A. John- 
son, Rhys. Rev. Lett. 62, 2825 (1989). 

[14] T. Raafat, J. P. Hulin, and H. J. Herrmann, Rhys. Rev. 
E 53, 4345 (1996). 

[15] K. L. Schick, and A. A. Verveen, Nature 251, 599 (1974). 

[16] X.-l. Wu, K. J. Maloy, A. Hansen, M. Ammi, and D. 
Bideau, Rhys. Rev. Lett. 71, 1363 (1993). 



[17] C. T. Vcje, and P. Dimon, Rhys. Rev. E 56, 4376 (1997). 
[18] P. Evcsquc and W. Mcftah, Int. J. Mod. Piiys. B 7, 1799 
(1993). 

[19] M. L. Hunt, R. C. Weathers, A. T. Lee, C. E. Brennen, 
and C. R. Wassgren, Rhys. Fluids 11 68 (1999). 

[20] C. R. Wassgren, M. L. Hunt, P. J. Freese, J. Palamara, 
and C. E. Brennen, Rhys. Fluids 14 3439 (2002). 

[21] K. Chen, M. B. Stone, R. Barry, M. Lohr, W. McConville, 
K. Klein, B. L. Sheu, A. J. Morss, T. Scheidemantel, and 
P. Schiffcr, Rhys. Rev. E 74, 011306 (2006). 

[22] H. Pacheco-Maxtinez, H.J. van Gerner, and J. C. Ruiz- 
Suarez, Rhys. Rev E 77 021303 (2008). 

[23] C. E. Davies and J. Foye, Trans. Inst. Chem. Eng. 69, 
269 (1991). 

[24] F. C. FrankUn and L. N. Johanson, Chem. Eng. Sci. 4, 

119 (1955). 

[25] C. D. Chitty and M. A. Spencer, Chemical Engineer- 
ing, Tripos Part 2. Research Project Report, University 
of Cambridge (1970). Nedderman informs us that this 
document has been discarded without regret (personal 
communication) . 

[26] C. S. Chang, H. H. Converse, and J. L. Steele, Trans. 
Am. Soc. Agr. Eng. 34, 1789 (1991). 

[27] J. F. Wambaugh, R. P. Behringer, J. V. Matthews, and 
P. A. Gremaud, Rhys. Rev. E 76, 051303 (2007) 

[28] U. Tiiziin, G. T. Houlsby, R. M. Nedderman, and S. B. 
Savage, Chem. Eng. Sci. 37, 1691 (1982). 

[29] J. Choi, A. KudroUi, and M. Z. Bazant, J. Phys.-Cond. 
Matt. 17, S2533 (2005). 

[30] C. S. O'Hern, L. E. Silbert, A. J. Liu, and S. R. Nagel, 
Rhys. Rev. E 68, 011306 (2003). 

[31] C. C. Thomas and D. J. Durian, now in progress. 



